<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">

<!-- 
	Copyright (C) 2007, 2008, 2009, 2010, 2011. PARP Research Group.
	<http://perception.inf.um.es>
	University of Murcia, Spain.

	This file is part of the QVision library.

	QVision is free software: you can redistribute it and/or modify
	it under the terms of the GNU Lesser General Public License as
	published by the Free Software Foundation, version 3 of the License.

	QVision is distributed in the hope that it will be useful,
	but WITHOUT ANY WARRANTY; without even the implied warranty of
	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
	GNU Lesser General Public License for more details.

	You should have received a copy of the GNU Lesser General Public
	License along with QVision. If not, see <http://www.gnu.org/licenses/>.
-->

<html><head><meta http-equiv="content-Type" content="text/html;charset=UTF-8">
<title>QVision: Qt&#39;s Image, Video and Computer Vision Library</title>
<meta name="title" content="QVision" />
<meta name="dc.title" content="QVision" />
<meta name="url" content="http://perception.inf.um.es/QVision" />
<meta name="author" content="PARP Research Group - http://perception.inf.um.es" />
<meta name="revisit-after" content="30 DAYS"/>
<meta name="robots" content="index,follow"/>
<meta name="classification" content="*">
<meta name="rating" content="Safe For Kids">
<meta name="distribution" content="GLOBAL"/>
<meta name="description" content="Qt's Image, Video and Computer Vision Library"/>
<meta name="page-topic" content="Computer Vision research and prototype programming"/>
<meta name="geo.country" content="ES" />

<!--
Keywords:
By license:		open source, gnu, lgpl, gpl, free
By theme:		computer vision, image processing, robotics, programming, source, development
By usage:		library, toolkit, framework, prototype, application
By programming specs:	object oriented, c++, block programming, reusability, gui, graphical, parallel computing, high performance, GPU, prototyping
Interoperability with:	Qt, GSL, GNU Scientific library, OpenCV, CGAL, QWT, CUDA, mplayer, IPP, Intel Image Performance Primitives, blas, lapack
Functionallity:		image features, matrix algebra, projective geometry, mser, function minimization, function optimization, canny operator, harris operator, corner detection, performance evaluation, cpu usage, graphical interface
Main data-types:	matrix, vector, tensor, quaternion, image, polyline
Video sources:		webcam, camera, stream
Devices:		embedded, desktop computer, laptop, mini-laptop
Authors:		PARP research group. University of Murcia, Spain.
-->

<meta name="keywords" content="augmented reality, sfm, structure from motion, open source, gnu, lgpl, gpl, free, computer vision, image processing, robotics, programming, source, development, library, toolkit, framework, prototype, application, object oriented, c++, block programming, reusability, gui, graphical, parallel computing, high performance, GPU, prototyping, Qt, GSL, GNU Scientific library, OpenCV, CGAL, QWT, CUDA, mplayer, IPP, Intel Image Performance Primitives, blas, lapack, image features, matrix algebra, projective geometry, mser, function minimization, function optimization, canny operator, harris operator, corner detection, performance evaluation, cpu usage, graphical interface, matrix, vector, tensor, quaternion, image, polyline, webcam, camera, stream, embedded, desktop computer, laptop, mini-laptop, University of Murcia, Spain, PARP research group, vision por computador"/>
<meta http-equiv="keywords" content="augmented reality, sfm, structure from motion, open source, gnu, lgpl, gpl, free, computer vision, image processing, robotics, programming, source, development, library, toolkit, framework, prototype, application, object oriented, c++, block programming, reusability, gui, graphical, parallel computing, high performance, GPU, prototyping, Qt, GSL, GNU Scientific library, OpenCV, CGAL, QWT, CUDA, mplayer, IPP, Intel Image Performance Primitives, blas, lapack, image features, matrix algebra, projective geometry, mser, function minimization, function optimization, canny operator, harris operator, corner detection, performance evaluation, cpu usage, graphical interface, matrix, vector, tensor, quaternion, image, polyline, webcam, camera, stream, embedded, desktop computer, laptop, mini-laptop, University of Murcia, Spain, PARP research group, vision por computador"/>
<meta http-equiv="pragma" content="no-cache"/>
<meta http-equiv="title" content="QVision"/>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
<link href="tabs.css" rel="stylesheet" type="text/css" />
<link rel="shortcut icon" href="favicon.ico" />
</head><body>

<table width="100%"><tr>
	<td><a href="http://perception.inf.um.es/"><img src="parp.png" border="0" /> <big>PARP Research Group</big></a></td>
	<td align="right"><a href="http://www.um.es/"><big>Universidad de Murcia</big> <img src="um.png" border="0" /></a></td>
</tr></table>

<hr /><br />

<table width="95%" align="center"><tr><td>

<!-- Generated by Doxygen 1.6.3 -->
<script type="text/javascript"><!--
var searchBox = new SearchBox("searchBox", "search",false,'Search');
--></script>
<div class="navigation" id="top">
  <div class="tabs">
    <ul>
      <li><a href="index.html"><span>Main&nbsp;Page</span></a></li>
      <li><a href="pages.html"><span>Related&nbsp;Pages</span></a></li>
      <li><a href="modules.html"><span>Modules</span></a></li>
      <li><a href="annotated.html"><span>Classes</span></a></li>
      <li class="current"><a href="files.html"><span>Files</span></a></li>
      <li>
        <div id="MSearchBox" class="MSearchBoxInactive">
          <form id="FSearchBox" action="search.php" method="get">
            <img id="MSearchSelect" src="search/search.png" alt=""/>
            <input type="text" id="MSearchField" name="query" value="Search" size="20" accesskey="S" 
                   onfocus="searchBox.OnSearchFieldFocus(true)" 
                   onblur="searchBox.OnSearchFieldFocus(false)"/>
          </form>
        </div>
      </li>
    </ul>
  </div>
  <div class="tabs">
    <ul>
      <li><a href="files.html"><span>File&nbsp;List</span></a></li>
      <li><a href="globals.html"><span>File&nbsp;Members</span></a></li>
    </ul>
  </div>
<h1>examples/OpenCV/delaunay/delaunay.cpp</h1><a href="delaunay_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*</span>
<a name="l00002"></a>00002 <span class="comment"> *      Copyright (C) 2007, 2008, 2009, 2010, 2011, 2012. PARP Research Group.</span>
<a name="l00003"></a>00003 <span class="comment"> *      &lt;http://perception.inf.um.es&gt;</span>
<a name="l00004"></a>00004 <span class="comment"> *      University of Murcia, Spain.</span>
<a name="l00005"></a>00005 <span class="comment"> *</span>
<a name="l00006"></a>00006 <span class="comment"> *      This file is part of the QVision library.</span>
<a name="l00007"></a>00007 <span class="comment"> *</span>
<a name="l00008"></a>00008 <span class="comment"> *      QVision is free software: you can redistribute it and/or modify</span>
<a name="l00009"></a>00009 <span class="comment"> *      it under the terms of the GNU Lesser General Public License as</span>
<a name="l00010"></a>00010 <span class="comment"> *      published by the Free Software Foundation, version 3 of the License.</span>
<a name="l00011"></a>00011 <span class="comment"> *</span>
<a name="l00012"></a>00012 <span class="comment"> *      QVision is distributed in the hope that it will be useful,</span>
<a name="l00013"></a>00013 <span class="comment"> *      but WITHOUT ANY WARRANTY; without even the implied warranty of</span>
<a name="l00014"></a>00014 <span class="comment"> *      MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the</span>
<a name="l00015"></a>00015 <span class="comment"> *      GNU Lesser General Public License for more details.</span>
<a name="l00016"></a>00016 <span class="comment"> *</span>
<a name="l00017"></a>00017 <span class="comment"> *      You should have received a copy of the GNU Lesser General Public</span>
<a name="l00018"></a>00018 <span class="comment"> *      License along with QVision. If not, see &lt;http://www.gnu.org/licenses/&gt;.</span>
<a name="l00019"></a>00019 <span class="comment"> */</span>
<a name="l00020"></a>00020 
<a name="l00045"></a>00045 <span class="preprocessor">#include &lt;stdio.h&gt;</span>
<a name="l00046"></a>00046 <span class="preprocessor">#include &lt;stdlib.h&gt;</span>
<a name="l00047"></a>00047 <span class="preprocessor">#include &lt;iostream&gt;</span>
<a name="l00048"></a>00048 <span class="preprocessor">#include &lt;math.h&gt;</span>
<a name="l00049"></a>00049 
<a name="l00050"></a>00050 <span class="preprocessor">#include &lt;QDebug&gt;</span>
<a name="l00051"></a>00051 <span class="preprocessor">#include &lt;QThread&gt;</span>
<a name="l00052"></a>00052 
<a name="l00053"></a>00053 <span class="preprocessor">#include &lt;QVApplication&gt;</span>
<a name="l00054"></a>00054 <span class="preprocessor">#include &lt;QVVideoReaderBlock&gt;</span>
<a name="l00055"></a>00055 <span class="preprocessor">#include &lt;QVDefaultGUI&gt;</span>
<a name="l00056"></a>00056 <span class="preprocessor">#include &lt;QVImageCanvas&gt;</span>
<a name="l00057"></a>00057 <span class="preprocessor">#include &lt;QVPolylineF&gt;</span>
<a name="l00058"></a>00058 
<a name="l00059"></a>00059 <span class="preprocessor">#include &lt;QVHarrisPointDetector&gt;</span>
<a name="l00060"></a>00060 <span class="preprocessor">#include &lt;qvip.h&gt;</span>
<a name="l00061"></a>00061 <span class="preprocessor">#include &lt;<a class="code" href="qvip_2qvimageio_8h.html" title="File from the QVision library.">qvimageio.h</a>&gt;</span>
<a name="l00062"></a>00062 
<a name="l00063"></a>00063 <span class="preprocessor">#ifndef DOXYGEN_IGNORE_THIS</span>
<a name="l00064"></a>00064 <span class="preprocessor"></span>
<a name="l00065"></a>00065 <span class="preprocessor">#ifdef OPENCV</span>
<a name="l00066"></a>00066 <span class="preprocessor"></span><span class="preprocessor">#include &lt;QVDelaunayTriangulation&gt;</span>
<a name="l00067"></a>00067 
<a name="l00068"></a>00068 QHash&lt;QPointF, QPointF&gt; openCVDelaunayTriangulation(<span class="keyword">const</span> QList&lt;QPointF&gt; &amp;selectedPoints)
<a name="l00069"></a>00069         {
<a name="l00070"></a>00070         <span class="comment">// STORAGE AND STRUCTURE FOR DELAUNAY SUBDIVISION</span>
<a name="l00071"></a>00071 
<a name="l00072"></a>00072         <span class="comment">// Set boundaries for the Delaunay triangulation algorithm.</span>
<a name="l00073"></a>00073         <span class="keywordtype">double</span>  minX = std::numeric_limits&lt;double&gt;::max(),
<a name="l00074"></a>00074                 minY = std::numeric_limits&lt;double&gt;::max(),
<a name="l00075"></a>00075                 maxX = std::numeric_limits&lt;double&gt;::min(),
<a name="l00076"></a>00076                 maxY = std::numeric_limits&lt;double&gt;::min();
<a name="l00077"></a>00077 
<a name="l00078"></a>00078         <span class="keywordflow">foreach</span>(QPointF point, selectedPoints)
<a name="l00079"></a>00079                 {
<a name="l00080"></a>00080                 minX = MIN(minX, point.x());
<a name="l00081"></a>00081                 minY = MIN(minY, point.y());
<a name="l00082"></a>00082                 maxX = MAX(maxX, point.x());
<a name="l00083"></a>00083                 maxY = MAX(maxY, point.y());
<a name="l00084"></a>00084                 }
<a name="l00085"></a>00085 
<a name="l00086"></a>00086         QVDelaunayTriangulation delaunayTriangulation(QRect(minX, minY, maxX - minX, maxY - minY));
<a name="l00087"></a>00087 
<a name="l00088"></a>00088         <span class="comment">// Add points to the Delaunay subdivision structure.</span>
<a name="l00089"></a>00089         <span class="keywordflow">foreach</span>(QPointF point, selectedPoints)
<a name="l00090"></a>00090                 delaunayTriangulation.addPoint(point);
<a name="l00091"></a>00091         <span class="keywordflow">return</span> delaunayTriangulation.getLinks();
<a name="l00092"></a>00092         }
<a name="l00093"></a>00093 
<a name="l00094"></a>00094 <span class="keyword">class </span>DelaunayBlock: <span class="keyword">public</span> <a class="code" href="classQVProcessingBlock.html" title="Base class for Block threads.">QVProcessingBlock</a>
<a name="l00095"></a>00095         {
<a name="l00096"></a>00096         <span class="keyword">private</span>:
<a name="l00097"></a>00097                 QList&lt;QVPolylineF&gt; map2polyline(<span class="keyword">const</span> QHash&lt;QPointF, QPointF&gt; &amp;map)
<a name="l00098"></a>00098                         {
<a name="l00099"></a>00099                         QList&lt;QVPolylineF&gt; polylines;
<a name="l00100"></a>00100                 
<a name="l00101"></a>00101                         QHashIterator&lt;QPointF, QPointF&gt; iterator(map);
<a name="l00102"></a>00102                 
<a name="l00103"></a>00103                         <span class="keywordflow">while</span> (iterator.hasNext())
<a name="l00104"></a>00104                                 {
<a name="l00105"></a>00105                                 iterator.next();
<a name="l00106"></a>00106                                 <a class="code" href="classQVPolylineF.html" title="Floating point polyline representation.">QVPolylineF</a> polyline;
<a name="l00107"></a>00107                                 polyline.append(iterator.key());
<a name="l00108"></a>00108                                 polyline.append(iterator.value());
<a name="l00109"></a>00109                                 polylines.append(polyline);
<a name="l00110"></a>00110                                 }
<a name="l00111"></a>00111                 
<a name="l00112"></a>00112                         <span class="keywordflow">return</span> polylines;
<a name="l00113"></a>00113                         }
<a name="l00114"></a>00114 
<a name="l00115"></a>00115         <span class="keyword">public</span>:
<a name="l00116"></a>00116                 DelaunayBlock(<span class="keyword">const</span> QString name): <a class="code" href="classQVProcessingBlock.html" title="Base class for Block threads.">QVProcessingBlock</a>(name)
<a name="l00117"></a>00117                         {
<a name="l00118"></a>00118                         addProperty&lt; QList &lt; QPointF &gt; &gt;(<span class="stringliteral">&quot;Feature locations&quot;</span>, inputFlag|outputFlag);
<a name="l00119"></a>00119                         addProperty&lt; QList&lt;QVPolylineF&gt; &gt;(<span class="stringliteral">&quot;Delaunay edges&quot;</span>, outputFlag);
<a name="l00120"></a>00120                         }
<a name="l00121"></a>00121 
<a name="l00122"></a>00122                 <span class="keywordtype">void</span> <a class="code" href="classQVProcessingBlock.html#a0283c799776e0c938538de0b0ae6b463" title="Function to be defined in subclasses with processing code.">iterate</a>()
<a name="l00123"></a>00123                         {
<a name="l00124"></a>00124                         <span class="comment">// 0. Read input property values.</span>
<a name="l00125"></a>00125                         <span class="keyword">const</span> QList&lt;QPointF&gt; selectedPoints = getPropertyValue&lt; QList&lt; QPointF &gt; &gt;(<span class="stringliteral">&quot;Feature locations&quot;</span>);
<a name="l00126"></a>00126                         <a class="code" href="classQVProcessingBlock.html#a61f390e8790a77175a6d6e1049a9d127" title="Stablish a performance breakpoint in the execution of the block.">timeFlag</a>(<span class="stringliteral">&quot;Read data&quot;</span>);
<a name="l00127"></a>00127 
<a name="l00128"></a>00128                         <span class="comment">// 1. Get Delaunay triangulation.</span>
<a name="l00129"></a>00129                         <span class="keyword">const</span> QHash&lt;QPointF, QPointF&gt; correspondences = openCVDelaunayTriangulation(selectedPoints);
<a name="l00130"></a>00130                         <a class="code" href="classQVProcessingBlock.html#a61f390e8790a77175a6d6e1049a9d127" title="Stablish a performance breakpoint in the execution of the block.">timeFlag</a>(<span class="stringliteral">&quot;Generated edges&quot;</span>);
<a name="l00131"></a>00131 
<a name="l00132"></a>00132                         <span class="comment">// 2. Generate the graph as a list of polylines for visualization.</span>
<a name="l00133"></a>00133                         <span class="keyword">const</span> QList&lt;QVPolylineF&gt; delaunayEdges = map2polyline(correspondences);
<a name="l00134"></a>00134                         <a class="code" href="classQVProcessingBlock.html#a61f390e8790a77175a6d6e1049a9d127" title="Stablish a performance breakpoint in the execution of the block.">timeFlag</a>(<span class="stringliteral">&quot;From correspondences to polylines&quot;</span>);
<a name="l00135"></a>00135 
<a name="l00136"></a>00136                         setPropertyValue&lt; QList&lt;QVPolylineF&gt; &gt;(<span class="stringliteral">&quot;Delaunay edges&quot;</span>, delaunayEdges);
<a name="l00137"></a>00137                         }
<a name="l00138"></a>00138         };
<a name="l00139"></a>00139 
<a name="l00140"></a>00140 <span class="preprocessor">#include &lt;QVYUV4MPEG2WriterBlock&gt;</span>
<a name="l00141"></a>00141 <span class="preprocessor">#include &lt;QVNumericPlot&gt;</span>
<a name="l00142"></a>00142 <span class="keywordtype">int</span> main(<span class="keywordtype">int</span> argc, <span class="keywordtype">char</span> *argv[])
<a name="l00143"></a>00143         {
<a name="l00144"></a>00144         <a class="code" href="classQVApplication.html" title="Main class for QVision applications.">QVApplication</a> app(argc, argv,
<a name="l00145"></a>00145                 <span class="stringliteral">&quot;Example program for QVision library. Obtains intrinsic and extrinsic camera parameters.&quot;</span>
<a name="l00146"></a>00146                 );
<a name="l00147"></a>00147 
<a name="l00148"></a>00148         <a class="code" href="classQVVideoReaderBlock.html" title="This is a generic video input processing block for the QVision architecture.">QVVideoReaderBlock</a> camera(<span class="stringliteral">&quot;Video&quot;</span>);
<a name="l00149"></a>00149 
<a name="l00150"></a>00150         <a class="code" href="classQVHarrisPointDetector.html">QVHarrisPointDetector</a> pointDetector(<span class="stringliteral">&quot;Harris corners detector&quot;</span>);
<a name="l00151"></a>00151         camera.linkProperty(&amp;pointDetector,<span class="stringliteral">&quot;Input image&quot;</span>);
<a name="l00152"></a>00152 
<a name="l00153"></a>00153         <span class="comment">// Camera calibrator.</span>
<a name="l00154"></a>00154         DelaunayBlock delaunayBlock(<span class="stringliteral">&quot;Delaunay&quot;</span>);
<a name="l00155"></a>00155         pointDetector.linkProperty(<span class="stringliteral">&quot;Feature locations&quot;</span>, delaunayBlock, QVProcessingBlock::SynchronousLink);
<a name="l00156"></a>00156 
<a name="l00157"></a>00157         <span class="comment">// Image outputs</span>
<a name="l00158"></a>00158         <a class="code" href="classQVImageCanvas.html" title="Display widget for QVImage objects and other data types.">QVImageCanvas</a> imageCanvas(<span class="stringliteral">&quot;Delaunay regions&quot;</span>);
<a name="l00159"></a>00159         pointDetector.linkProperty(<span class="stringliteral">&quot;Input image&quot;</span>, imageCanvas);
<a name="l00160"></a>00160         delaunayBlock.linkProperty(<span class="stringliteral">&quot;Feature locations&quot;</span>, imageCanvas);
<a name="l00161"></a>00161         delaunayBlock.linkProperty(<span class="stringliteral">&quot;Delaunay edges&quot;</span>, imageCanvas);
<a name="l00162"></a>00162         imageCanvas.setColor(<span class="stringliteral">&quot;Delaunay edges&quot;</span>, Qt::red);
<a name="l00163"></a>00163 
<a name="l00164"></a>00164         <a class="code" href="classQVDefaultGUI.html" title="Main user interface widget for QVision applications.">QVDefaultGUI</a> gui;
<a name="l00165"></a>00165 
<a name="l00166"></a>00166         <span class="keywordflow">return</span> app.exec();
<a name="l00167"></a>00167         }
<a name="l00168"></a>00168 <span class="preprocessor">#else</span>
<a name="l00169"></a>00169 <span class="preprocessor"></span><span class="keywordtype">int</span> main(<span class="keywordtype">int</span> argc, <span class="keywordtype">char</span> *argv[])
<a name="l00170"></a>00170         {
<a name="l00171"></a>00171         std::cout &lt;&lt; <span class="stringliteral">&quot;ERROR: OpenCV compatibility was not activated in QVision compilation.&quot;</span> &lt;&lt; std::endl;
<a name="l00172"></a>00172         <span class="keywordflow">return</span> -1;
<a name="l00173"></a>00173         }
<a name="l00174"></a>00174 <span class="preprocessor">#endif  // OPENCV</span>
<a name="l00175"></a>00175 <span class="preprocessor"></span><span class="preprocessor">#endif  // DOXIGEN_IGNORE_THIS</span>
<a name="l00176"></a>00176 <span class="preprocessor"></span>
</pre></div></div>
</td></tr></table>

<br /><hr><br />
<center><a href="http://perception.inf.um.es/QVision">QVision framework</a>.
<a href="http://perception.inf.um.es">PARP research group</a>.
Copyright &copy; 2007, 2008, 2009, 2010, 2011.</center>
<br />
</body>
</html>
